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We extend the twistor string theory inspired formalism introduced in hep-th/0407214 for 
calculating loop amplitudes in M = 4 super Yang-Mills theory to the case of M = 1 (and 
M =2) super Yang-Mills. Our approach yields a novel representation of the gauge theory 
amplitudes as dispersion integrals, which are surprisingly simple to evaluate. As an ap- 
Ph! plication we calculate one-loop maximally helicity violating (MHV) scattering amplitudes 
with an arbitrary number of external legs. The result we obtain agrees precisely with the 
expressions for the M = 1 MHV amplitudes derived previously by Bern, Dixon, Dunbar 



and Kosower using the cut-constructibility approach. 
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1 Introduction 



In a remarkable paper [1] it was proposed that perturbative J\f = 4 super Yang-Mills 
(SYM) is dual to the topological B model on super twistor space CP 3 ' 4 . Interestingly, 
this duality relates the perturbative expansion of gauge theory amplitudes to a Dl-brane 
instanton expansion on the string theory side. The relevant instantons correspond to 
algebraic curves embedded holomorphically in super twistor space. Their degree d and 
genus g are related to the number of negative helicity gluons q and number of loops I of 
the amplitude as d = q — 1 + I and g < I. 

In [1] the maximally helicity violating (MHV) tree amplitudes [2,3] were reproduced 
directly from a computation in the B model with d — 1, g — curves. For more gen- 
eral amplitudes the story becomes more involved, since in [1] it was already pointed out 
that in principle both connected and disconnected instantons can contribute to a given 
amplitude. The approach using connected instantons was pursued further in [4-6], and 
agreement with existing results was found. On the other hand, [7] introduced a new dia- 
grammatic method to calculate tree-level amplitudes using MHV amplitudes as effective 
vertices after a suitable off-shell continuation. This method is related to the string the- 
ory approach using contributions from completely disconnected instantons only, and is 
extremely efficient for calculating tree amplitudes. The authors of [8] argued that both 
computations are equivalent, and that the instanton contributions localize on singular 
curves corresponding to intersecting "lines" , i.e. degree one curves in twistor space. 

The method of [7] led to a rederivation of previously known results in a much faster 
way, and also allowed the calculation of new scattering amplitudes with increasing helicity 
violation as well as with fermionic and scalar external particles [9]-[16]. The next logical 
question to ask was clearly whether this procedure could be extended to the calculation 
of amplitudes at loop level. This question was positively answered by three of the present 
authors in [17], where it was shown how to combine MHV vertices into loop diagrams in 
Af = 4 SYM theory. The result was a new representation of the scattering amplitudes 
in terms of dispersive integrals which, rather surprisingly, proved to be tractable, and led 
to perfect agreement with the expressions for the one-loop MHV scattering amplitudes 
in M = 4 SYM previously obtained by Bern, Dixon, Dunbar and Kosower (BDDK) 
in [18] using the cut-constructibility approach. As a bonus, the analysis of [17] gave a 
novel representation of the "easy two-mass" (2me) scalar box function, containing one 
less dilogarithm and one less logarithm than the traditional one of BDDK. 

The twistor string theory originally proposed in [1] and other versions proposed in 
[19,20] cannot reproduce gauge theory amplitudes at loop level [21], due to the fact 
that conformal supergravity fields propagate in the loops. Nevertheless, one can study 
known loop amplitudes to extract information about their localisation properties in twistor 
space. This was done in [22] for supersymmetric MHV one-loop amplitudes (and some 
non-supersymmetric amplitudes), and localisation was indeed found onto three types of 
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diagrams. The twistor space picture suggested by the result of [17] differed from that 
of [22], in that one of the class of diagrams proposed in the analysis of [22] is absent. 
The resolution of this discrepancy was recently [23] linked to a subtlety in the use of the 
differential operators employed to establish localisation, and led to the appearance of a 
"holomorphic anomaly". With this taken into account, it was shown in [23] that, as far 
as the one-loop Af = 4 SYM amplitudes are concerned, the expectations from [17] are 
confirmed and the one-loop MHV amplitudes in Af = 4 localise on pairs of lines in twistor 
space which are joined by two twistor space propagators. The holomorphic anomaly was 
evaluated explicitly in [24,25], and twistor space localisation was further elucidated in [24] 
by transforming the results of [7] and [17] directly to twistor space. Furthermore, the holo- 
morphic anomaly was recently exploited as a new tool for deriving one-loop next-to-MHV 
amplitudes in Af = 4 SYM [25, 26] in combination with the powerful cut-constructibility 
approach. These and other new one-loop amplitudes in Af = 4 SYM were recently derived 
in [27] using cut-constructibility. 

It now appears plausible that the entire quantum theory of Af = 4 SYM will have a 
description based upon the MHV diagram approach, which in turn reflects the properties 
of localisation in twistor space. Moreover, we now have direct methods to test this, 
using MHV vertices assembled into MHV diagrams according to well-defined rules. We 
expect further evidence to emerge with the continuing study of non-MHV tree and one- 
loop amplitudes, as well as higher-loop amplitudes. It should be fairly straightforward to 
check that the few remaining existing results that have so far not been obtained by these 
new methods are reproduced. 

We find it intriguing that, although the twistor string/gauge theory correspondence 
is spoiled at loop level, one can still use MHV vertices and combine them to obtain the 
correct Af = 4 SYM loop amplitudes [17]. Perhaps even more surprisingly, we will show 
in this paper that the applicability of the method of [17] actually holds for theories with 
less supersymmetry as well, which are interesting not least for potential applications to 
phenomenology. Specifically, in this paper we will obtain one-loop MHV amplitudes in 
Af = 1 SYM by combining MHV vertices into one-loop diagrams following the procedure 
proposed in [17]. The result we find agrees perfectly with the BDDK computation of [28]. 2 

The rest of this paper is organised as follows. In Section 2 we begin by reviewing the 
expression for the one-loop MHV scattering amplitudes in Af = 1 SYM derived by BDDK 
in [28]. There, we also present a slightly simplified version of the BDDK result which 
will be useful in making the comparison with the results derived from MHV vertices. 
Section 3 reviews relevant aspects of the work of [17], thereby establishing the set-up 
for the calculation of loop amplitudes with MHV vertices. In Section 4 we turn to the 
formulation and explicit calculation of the one-loop MHV amplitudes in Af = 1 SYM. 
Finally, in three Appendices we summarise some technical results which are required in 
the course of the calculations presented in Section 4. 

2 The arguments presented here also apply directly to TV = 2 SYM. 
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2 The J\f = 1 MHV amplitudes at one loop 



The expression for the MHV amplitudes at one loop in M = 1 SYM was obtained for 
the first time by BDDK in [28], using the cut-constructibility method. We will shortly 
give their explicit result, and then simplify it by introducing appropriate functions. This 
turns out to be useful when we compare the BDDK result to that which we will derive 
by using MHV diagrams. 

In order to obtain the one-loop MHV amplitudes in M = 1 and M = 2 SYM it 
is sufficient to compute the contribution ^4^ =1,chiral to the one-loop MHV amplitudes 
coming from a single M = 1 chiral multiplet. This was calculated in [28], and the result 
turns out to be proportional to the Parke- Taylor MHV tree amplitude [2] 
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as is also the case with the one-loop MHV amplitudes in M = 4 SYM. However, in 
contradistinction with that case, the remaining part of the M = 1 amplitudes depends 
non-trivially on the position of the negative helicity gluons % and j. The result obtained 
in [28] is: 
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where if := (pi + pi + \ + ■ ■ ■ + p i+ k-i) 2 for k > 0, and tf^ = t'f for k < 0. The sums in 
the second line of ()2.2j) cover the ranges C m and V m defined by 
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+ 1, 1}, 
[ {j + 1,j + 2,...,*-1}, 
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where g rjS := J^ s =r fcj. Notice that both coefficients b l ^ s and c^ a are symmetric under the 
exchange of i and j. In the case of b this is evident; for c, this is also manifest as it is 
expressed as the product of two antisymmetric quantities. The function B in the first line 
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Figure 1: The box function F of (12.7)1 . whose finite part B, Eq. (|2.9)l . appears in the 
M = 1 amplitude ()2.2jl . T/ie iu>o external gluons with negative helicity are labelled by i 
and j . The legs labelled by s and m correspond to the null momenta p and q respectively in 
the notation of ([2.9)1 . Moreover, the quantities tj^+i , tm~ m ^ tm+i 1 " 1 ^ ^I+i S appearing 
in the box function B in (|2.19|) correspond to the kinematical invariants t := (Q +p) 2 , 
s : = (P + p) 2 ; Q 2 ; p 2 in the notation of (|23)l . witt p + g + P + Q = 0. 

of ()2.2j) is the "finite" part of the easy two-mass (2me) scalar box function F(s, t, P 2 , Q 2 ), 
with 

1 



F(s,t : P 2 ,Q 2 ' 



(-*)- + (-t)--(-P^--(-Q^- + P( M ,P 2 ,Q 2 ). (2.7) 
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Figure 2: A triangle function, corresponding to the first term T e (p m , g a+ i jm _i, q m +i, a ) in 
the second line of ()2.19|) . p, Q and P correspond to p m , q m +i,a and g a+ i j?n _i m t/ie notation 
of Eq. (|2.19|) . w;/iere j e Q, i <E P. In particular, Q 2 — > tj^.i and P 2 — > tm~ m+1 ^ . 



As in [17], we have introduced the following convenient kinematical invariants: 

8 := (P + p) 2 , t := (P + g) 2 , (2.8) 

where p and g are null momenta, and P and Q are in general massive, with p + q + P + Q = 
0. 3 In [17] the following new expression for B was found: 

B(s, t, P 2 , Q 2 ) = Li 2 (l - aP 2 ) + Li 2 (l - aQ 2 ) - Li 2 (l - as) - Li 2 (l - at) , (2.9) 

where 

P 2 + Q 2 -s-t . 

The expression ()2.9)1 contains one less dilogarithm and one less logarithm than the tradi- 
tional form used by BDDK, 

P 2 \ / P 2 \ / Q 2 \ ( Q 2 

r) +Li2 ( 1 -T) +Li2 ( 1 -i 



P(s,t,P 2 ,Q 2 ) = Li 2 (l - — J + Li 2 (l- V)+Li 2 (l-— ) + Li 2 (l-^- 

P 2 Q 2 \ , 1, 2(S 



Li2(l -^) + iv(f). (,n; 



3 The kinematical invariant s = (P +p) 2 should not be confused with the label s which is also used to 
label an external leg (as in Figure 1 for example). The correct meaning will be clear from the context. 
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Figure 3: This triangle function corresponds to the second term in the second line of 1)2.19)1 
- where i and j are swapped. As in Figure 2, p, Q and P correspond to p m , q m +i, a and 
Qa+i,m-i i n the notation of Eq. (12.19)1 . where nowi G Q, j G P. In particular, Q 2 — > 

and p2 ^ t [m-a-l] 

The agreement of ()2.9)1 with ()2.11jl was discussed and proved in Section 5 of [17]. 4 In 
Figure 1 we give a pictorial representation of the box function F defined in (J2.7j) (with 
the leg labels identified by s — > p, m — > q). 

Finally, infrared divergences are contained in the bubble functions Ko(t), defined by 

: = ' (2 - 12) 

We notice that in order to re-express (|2.2|) in a simpler form, it is useful to introduce the 
triangle function [22] 

T(p,P,Q) := Q2 _ p2 , (2-13) 
with p + P + Q = 0. A diagrammatic representation of this function is given in Figure 2 



4 More precisely, this agreement holds only in certain kinematical regimes e.g. in the Euclidean region 
where all kinematical invariants are negative. More care is needed when analytically continuing the 
amplitude to the physical region. The usual prescription of replacing a kinematical invariant s by s + ie 
and continuing s from negative to positive values gives the correct result only for our form of the box 
function 1)2.9(1 . whereas (|2.11(l has to be amended by correction terms [29]. 
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(with m + — > p). We also find it useful to introduce an e-dependent triangle function 
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Te( P ,P,Q) ■= ~ e { _ ^P' 6 ■ (2.14) 

As long as P 2 and Q 2 are non-vanishing, one has 

UmT e (p,P,Q) = T(p,P,Q), P V , Q V . (2.15) 

If either of the invariants vanishes, one has a different limit. For example, if Q 2 = 0, one 
has 2 

T e (p,P,Q)\ Q2=0 — » -; ( ~p 2 } \ e^O. (2.16) 
We will call these cases "degenerate triangles". 

The usefulness of the previous remark stems from the fact that precisely the quantity 
(1/e)- [(— P 2 )~ e /P 2 ] appears in the last line of ()2.2|) - the bubble contributions. Therefore, 
these can be equivalently obtained as degenerate triangles, i.e. triangles where one of the 
massive legs becomes massless. 

Specifically, we notice that the four degenerate triangles (bubbles) in the last line of 
(|2.2j) can be precisely obtained by including the "missing" index assignments in T> m and 
C ■ 

(m = i + 1, a = i — 1) , (m = j — 1, a = j) for V m , (2-17) 
which correspond to two degenerate triangles, and 

(m — j + 1, a — j — 1) , (m = i — 1, a — i) for C m , (2-18) 

corresponding to two more degenerate triangles. 

In conclusion, the previous remarks allow us to rewrite ()2.2|) in a more compact form 
as follows: 

j-i i-l 



jAT=l,chiral _ jtree I uj D('+[ s " m ] +[s-m] .[s-m-1] ,[m-s-l]\ /n i q\ 

•^n — ' \ / j / j u m,s JD \ L m+l i l m i h m+\ i b s+l I \^- Ly ! 

^m=i+l s=7+l 

XI C m,a^(Pm>9a+l,m-l,gm+l,a) + (« < > j) >■ 



1 

+ 



1 - 2e 



m=i+l a=j 



In the previous expression it is understood that we only keep terms that survive in the 
limit e — > 0. This means that the factor 1/(1 — 2e) can be replaced by 1 whenever the term 
in the sum is finite, i.e. whenever the triangle is non-degenerate. However, in the case of 



5 The function T e (p, P, Q) denned in l|2.14|l arises naturally in the twistor- inspired approach which will 
be developed in Sections 3 and 4. 
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degenerate triangles, which contain infrared divergent terms, we have to expand this factor 
to linear order in e. In the notation of ([2. 19)1 . q^+i a = an d Qa+i m-i = tm~ m+1 ^; in 

Figure 2, these invariants correspond to Q 2 and P 2 respectively, where j G Q, i G P. In 
the sum with % <-> j, one would have g^ +la = v^ t _ 1 , m-i = ^i+i a , corresponding 
respectively to Q 2 and P 2 in Figure 3, with i G Q, j G P. It is the expression (j2.19j) for 
the Af = 1 chiral multiplet amplitude which we will derive using MHV diagrams. 



3 A brief review of loop diagrams from MHV vertices 

In this Section we briefly review the method proposed in [17] to compute loop amplitudes 
in super symmetric gauge theories, referring the reader to Sections 3-5 of that paper for 
a more detailed discussion of this procedure, as well as for an example of its application 
to the calculation of one-loop MHV amplitudes in M = 4 SYM. 

Before starting, let us first recall that in what follows we will be dealing with the so- 
called colour-stripped or partial amplitudes. More specifically, for an n-particle scattering 
amplitude we will compute the term proportional to Tr(T ai ■ • -T an ), where the T a 's are 
the generators of the gauge group. The full planar amplitude is then obtained by summing 
over non-cyclic permutations. It is a remarkable result of BDDK that, at one loop, non- 
planar amplitudes are simply obtained as a sum over permutations of the planar ones. 
This is discussed in Section 7 of [18], where it was also noted that this applies to a generic 
theory with adjoint particles running in the loops, such as Af = 1, 2, 4 SYM. At the level 
of group theory factors, the diagrammatics for building one-loop MHV amplitudes in the 
BDDK cut-constructibility approach and in the approach discussed in this paper (as well 
as in [17], for the Af = 4 case) are precisely the same. Hence, the agreement at one loop 
between the two methods at planar level in fact trivially extends to subleading corrections 
in 1/JV. 

We now come back to the description of the approach of [17] to loop amplitudes. The 
procedure used there was: 

1. Lift the MHV tree- level scattering amplitudes to vertices, by continuing the internal 
lines off shell using a prescription equivalent to that of CSW. Internal lines are then 
connected by scalar propagators which join particles of the same spin but opposite 
helicity. 

2. Build MHV diagrams with the required external particles at loop level using the 
MHV tree-level vertices, and sum over all independent diagrams obtained in this 
fashion for a fixed ordering of external helicity states. 

3. Re-express the loop integration measure in terms of the off-shell parametrisation 
employed for the loop momenta. 
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4. Analytically continue to 4 — 2e dimensions in order to deal with infrared divergences, 
and perform all loop integrations. 



Using this method, we will show in Section 4 that combining MHV vertices into one-loop 
diagrams precisely yields the results for the contribution of the chiral multiplet to MHV 
amplitudes, as obtained by BDDK using the cut-constructibility approach. 

We start off by discussing the off-shell continuation which was used in [17], and found 
to be very useful for calculating loop diagrams. Consider a generic off-shell momentum 
vector, L. On general grounds, it can be decomposed as [13, 14] 

L = I + zt] , (3.1) 

where I 2 = 0, and 77 is a fixed and arbitrary null vector, rf = 0; z is a real number. 
Equation (|3.1J) determines Z clS db function of L: 



L 2 



2(L V ) ■ 



(3.2) 



Using spinor notation, we can write I and r\ as l a a = l a l a , r/ a a = r\ a r}u- It then follows 
that 6 

T ri« 

L = ^r, (3.3) 



n a T ■ 

l a . [6 A) 

We notice that ()3.3|) and ()3.4j) coincide with the CSW prescription proposed in [7] to 
determine the spinor variables / and / associated with the non-null, off-shell four-vector 
L defined in (13. 1|) . The denominators on the right hand sides of (j3.3j) and ()3.4j) turn out 
to be irrelevant for our applications, since the expressions we will be dealing with are 
homogeneous in the spinor variables l a ; hence we will discard them. 

To proceed further, we need to re-express the usual integration measure d A L over the 
loop momentum L in terms of the new variables / and z introduced previously. We found 
that 7 

d 4 L dj\f(l) dz 

L 2 + ie ~ Ai z + is ' 
where we have introduced the Nair measure [31] 



dtf(l) := (I dl)d 2 l - [I dl] d 2 l . (3.6) 



6 Spinor inner products are defined as (A/i) := e a ij\ a ^^ [A p\ := e^A"/!^. 

7 The ie prescription in the left and right hand sides of 1)3. 5|l was understood in [17], and, as stressed 
in [f3,30], it is essential in order to correctly perform loop integrations. 



9 



Eq. (J3.5j) is key to the procedure. It is important to notice that the product of the measure 
factor with a scalar propagator d A L/ (L 2 +ie) in ()3.5|) is independent of the reference vector 
rj. In [31], it was noticed that the Lorentz invariant phase space measure for a massless 
particle can be expressed precisely in terms of the Nair measure: 

d H5^\l 2 ) = , (3.7) 

where, as before, we write the null vector I as l ad , = l a la, and in Minkowski space we 
identify 1 = 1*. 



Next, we observe that in computing one-loop diagrams, the four-dimensional integra- 
tion measure which appears is 

dM := ^^5^(L 2 -L 1 + P L ), (3.8) 
L x L 2 

where L x and L 2 are loop momenta, and Pl is the external momentum flowing outside 
the loop 8 so that L 2 — L\ + Pl — 0. Next we express L x and L 2 as in (|3.1|) . 

Li- t a,a ^ialia ~\~ %i TjaVot 5 ^ 1,2. (3.9) 

Using ([3.9)1 . we rewrite the argument of the delta function as 

L 2 -L 1 + P L = l 2 -l 1 + P L , z , (3.10) 

where we have defined 

P L;z ■■= Pl - zv , (3.H) 

and 

z := Z\ — z 2 . (3-12) 

Notice that we use the same r] for both the momenta Li and L 2 . Using ([3.9)1 . we can 
finally recast ()3.8j) as [17] 

dM = dzidz2 dup _^ 
Z\ z 2 

where 

dUPS(l 2 ,-h-P L]Z ) := d%5^\ll) d%5 {+ \ll) 5 {4 \l 2 -h + P L ; Z ) (3.14) 

is the two-particle Lorentz invariant phase space (LIPS) measure. The integration mea- 
sure dM. as it is expressed on the right hand side of ([3.13)) can now be immediately 
dimensionally regularised; this is accomplished by simply replacing the four-dimensional 
LIPS measure by its continuation to D = 4 — 2e dimensions, 

d D UPS(l 2 ,-h;P Llz ) := d%5M(l 2 1 ) d%5^(l 2 2 ) 5W(l 2 -h+P L ., g ) . (3.15) 



3 In our conventions, all external momenta are outgoing. 
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Eq. ()3.13j) was one of the key results of [17]. It gives a decomposition of the original 
integration measure into a .D-dimensional phase space measure and a dispersive measure. 
According to Cutkosky's cutting rules [32], the LIPS measure computes the discontinuity 
of a Feynman diagram across its branch cuts. Which discontinuity is evaluated is deter- 
mined by the argument of the delta function appearing in the LIPS measure; in (13.13)1 this 
is Pl,z (defined in (J3.ll)) ). Notice that Pl- z always contains a term proportional to the 
reference vector r], as prescribed by ()3.11j) . Finally, discontinuities are integrated using 
the dispersive measure in (|3.13|) . thereby reconstructing the full amplitude. 

As a last remark, notice that, in contradistinction with the cut-constructibility ap- 
proach of BDDK, here we sum over all the cuts - each of which is integrated with the 
appropriate dispersive measure. 

4 MHV one-loop amplitudes in Af = 1 SYM from 
MHV vertices 

In the last section we reviewed how MHV vertices can be sewn together into one-loop 
diagrams, and how a particular decomposition of the loop momentum measure leads to a 
representation of the amplitudes strikingly similar to traditional dispersion formulae. This 
method was tested successfully in [17] for the case of MHV one-loop amplitudes in M = 4 
SYM. In the following we will apply the same philosophy to amplitudes in M = 1 SYM, 
in particular to the infinite sequence of MHV one-loop amplitudes, which was obtained 
using the cut-constructibility approach [28], and whose twistor space picture has more 
recently been analysed in [22]. 

Similarly to the Af = 4 case, the one-loop amplitude has an overall factor proportional 
to the MHV tree-level amplitude, but, as opposed to the Af = 4 case, the remaining 
one-loop factor depends non-trivially on the positions % and j of the two external negative 
helicity gluons. This is due to the fact that a different set of fields is allowed to propagate 
in the loops. 

The MHV Feynman diagrams contributing to MHV one-loop amplitudes consist of 
two MHV vertices connected by two off-shell scalar propagators. If both negative helicity 
gluons are on one MHV vertex, only gluons of a particular helicity can propagate in the 
loop. This is independent of the number of supersymmetries. On the other hand, for 
diagrams with one negative helicity gluon on one MHV vertex and the other negative 
helicity gluon on the other MHV vertex, all components of the supersymmetric multi- 
plet propagate in the loop. In the case of J\f = 4 SYM this corresponds to helicities 
h = —1, —1/2,0, 1/2, 1 with multiplicities 1,4,6,4, 1, respectively; for the Af = 1 vector 
multiplet the multiplicities are 1, 1, 0, 1, 1. Hence, we can obtain the Af = 1 amplitude by 
simply taking the Af = 4 amplitude and subtracting three times the contribution of an 
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Af = 1 chiral multiplet, which has multiplicities 0, 1, 2, 1, 0. 9 



This supersymmetric decomposition of general one-loop amplitudes is useful as it 
splits the calculation into pieces of increasing difficulty, and allows one to reduce a one- 
loop diagram with gluons circulating in the loop to a combination of an Af = 4 vector 
amplitude, an Af = 1 chiral amplitude and finally a non-supersymmetric amplitude with 
a scalar field running in the loop. 

In our case, the supersymmetric decomposition takes the form 

^Af=l,vector j^N=A g ^A/"=l, chiral j\ 

where n denotes the number of external lines. Since the Af = 4 contribution is known, one 
needs to determine ^4.^ =1 > chiral using MHV diagrams. To be precise, we are only addressing 
the computation of the planar part of the amplitudes. However this is sufficient, since 
at one-loop level the non-planar partial amplitudes are obtained as appropriate sums 
of permutations of the planar partial amplitudes [18], as discussed at the beginning of 
Section 3. 




Figure 4: A one-loop MHV diagram, computed in 4^4 using MHV amplitudes as interac- 
tion vertices, with the CSW off-shell prescription. The two external gluons with negative 
helicity are labelled by i and j. 



Therefore our task consists of 

1. Evaluating the class of diagrams where we allow all the helicity states of a chiral 
multiplet, 

he {-1/2, 0,0, 1/2}, (4.2) 
to run in the loop. We depict the prototype of such diagrams in Figure 4. 



5 We can also obtain the M = 2 amplitude in a completely similar way. 
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2. Summing over all diagrams such that each of the two MHV vertices always has one 
external gluon of negative helicity. Assigning i~ to the left and j~ to the right, the 
summation range of mi and 777.2 is determined to be: 

j + 1 < 777! < i , i < 777 2 < j — 1 . (4.3) 

Therefore we get 

^A/=l,chiral = JdMA(-h, mi ,...,r,...,m 2 ,l 2 ) 

roi,ms,/i 

■A(-h,m 2 + l,...,j~,...,m 1 -l,l 1 ) , (4.4) 

where the summation ranges of h, mi and m 2 are given in ()4.2j) . ()4.3j) . Notice that, in 
order to compute the loop amplitude ()4.4j) . we make use of the integration measure dAi, 
given in (j3.13|) . which was found in [17]. 

After some spinor algebra and after performing the sum over the helicities h, the 
integrand of f)4.4j) becomes 

_ • .tree (m 2 (m 2 +l)) ((7771 -1) 777!) (j lj) (j h) {i h) {j k) ( . 

n ' (^') 2 (^i^i)(K-l)^i)(^2/ 2 )((m 2 + l)/ 2 ) • { - } 

The focus of the remainder of this section will be to evaluate the integral in f)4.4j) explicitly. 
Since — factors out completely, we will drop it and only reinstate it at the very end 
of the calculation. 

The integrand (without this factor) can be rewritten in terms of dot products of 
momentum vectors, 

J\f 

1 = (7 • jf (m! • h) (K-l) • k) (m 2 • l 2 ) ((m 2 + l) • l 2 ) ' (46) 

with 

U = tr + (f ^mx-l^m xhMi) tr + (hK ztymz+lhMi) ■ ( 4 ' 7 ) 

Af is a product of Dirac traces, where the tr + symbol indicates that the projector (l+7 5 )/2 
has been inserted. 

Next, notice that each of these Dirac traces involving six momenta can be expressed 
in terms of simpler Dirac traces involving only four momenta. For the first factor of (|4.7|) 
we find 

tr + (J f 1 # mi _ 1 # fm J f 1 jyfi) = 2(mi •Zi)tr + (^ i ^ mi _ifi) - 2((m 1 -l) • Z x )tr + (^-^/i) , 

(4.8) 

where 

tr+ {tyaM4d) = 2 [(a ■ 6) (c • d) - (a ■ c) (6 ■ d) + (a • d) (6 ■ c)] - 2ze(a, 6, c, d) . (4.9) 
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The second factor in (|4.7jl takes a similar form. Consequently, the integrand becomes a 
sum of four terms, one of which is 

tr + ($iPjjpmJ_ l) tr + (jCjjCjjj/mgf 2) ('Ain'l 

(z-i) 2 K^i)K-z 2 ) ' 1 ' 1 

The other three terms are obtained by replacing m\ with mi — 1 and/or 777.2 with 7712 + 1 
in (j4.1()j) and come with alternating signs. Note that the original expression (|4.5jl is 
symmetric in i, and j, although when we make use of the decomposition (|4.10|) this 
symmetry is no longer manifest. We will symmetrize over i and j at the end of the 
calculation in order to make this exchange symmetry manifest in the final expression. 

In the next step we have to perform the phase space integration, which is equivalent to 
the calculation of a unitarity cut with momentum P Lz = Y^ mi ^i~ zr l flowing through the 
cut. Note that, as explained in Section 3, the momentum is shifted by a term proportional 
to the reference momentum 77. The term (l\ ■ mi)(l 2 • rn 2 ) i n the denominator in ()4.10|) 
corresponds to two propagators, hence the denominator by itself corresponds to a cut box 
diagram. However, the numerator of (j4.1()j) depends non-trivially on the loop momentum, 
so that in fact ()4.10|) corresponds to a tensor box diagram, not simply a scalar box diagram. 
Using the Passarino-Veltman method [33], we can reduce the expression (j4.10J) . integrated 
with the LIPS measure, to a sum of cuts of scalar box diagrams, scalar and vector triangle 
diagrams, and scalar bubble diagrams. This procedure is somewhat technical, and details 
are collected in Appendix A. Luckily, the final result takes a less intimidating form than 
the intermediate expressions. We will now present the result of these calculations after 
the LIPS integration. 

We first observe that loop integrations are performed in 4 — 2e dimensions. It turns 
out that singular 1/e terms appearing at intermediate steps of the phase space integration 
cancel out completely. Notice that this does not mean that the final result will be free 
of infrared divergences. In fact the dispersion integral can and does give rise to 1/e 
divergent terms but there cannot be any 1/e 2 terms, as expected for the contribution of 
a chiral multiplet [28] . The 1/e divergences in the scattering amplitude correspond to the 
bubble contributions in ()2.2|) . or degenerate triangles contributions in ()2.19j) . as explained 
in Section 2. In Appendix A we show that the finite terms of the phase space integral 
combine into the following simple expression: 

C = C{m\ — l,m 2 ) — C(mi,m 2 ) + C(mi,m 2 + 1) — C(m\ — l,m 2 + 1) , (4.11) 
with 10 



C(mi, m 2 ) 



27T (Pi ^ 



L;zJ 



1 — 2e (i ■ j) 2 (mi • 7772) 
2irT(m l , m 2 ,m 2 



T(m l ,m 2 ,P L ;z) T(m 2 ,m 1 ,P L . z 



{mi ■ P L . Z ) (m 2 ■ P L . Z 



(i- j) 2 {mi -m 2 ) 2 



(P 2 J- e log(l-a z P 2 2 ) , (4.12) 



10 



In 14.12J1 we omit an overall, finite numerical factor that depends on e. This factor, which can be 



read off from (IB . 1 21 . is irrelevant for our discussion. 



14 



where 



T(m 1 ,m 2 , P) 



tl_|_ ^i^j^ mi ^ tr+ (|»i^j^m2^mi) j 



mi • m 2 



(4.13) 



a 



z 



and 



N(P) := (mi ■ m 2 )P 2 - 2(mi ■ P)(m 2 ■ P) . 



(4.14) 



A closer inspection of 1)4.12)1 reveals that the first line of that expression corresponds to two 
cuts of scalar triangle integrals, up to an e-dependent factor and the explicit ^-dependence 
of the two numerators. The second line is a term familiar from [17], corresponding to the 
P|. z -cut of the finite part B of a scalar box function, defined in (J2.9j) (see also (|2.7j) ). 
The full result for the one-loop MHV amplitudes is obtained by summing over all possible 
MHV diagrams, as specified in (|Q|) and ([OJl . ([Of . 

We begin our analysis by focusing on the box function contributions in ()4.12|) . and 
notice the following important facts: 

1. By taking into account the four terms in ()4.11|) and summing over Feynman dia- 
grams, we see each fixed finite box function B appears in exactly four phase space 
integrals, one for each of its possible cuts, in complete similarity with [17]. It was 
shown in Section 5 of that paper that the corresponding dispersion integration over 
z will then yield the finite B part of the scalar box functions P. It was also noted 
in [17] that one can make a particular gauge choice for r\ such that the ^-dependence 
in iV disappears. This happens when rj is chosen to be equal to one of the massless 
external legs of the box function. The question of gauge invariance is discussed in 
Appendix C. 

2. The coefficient multiplying the finite box function is precisely equal to 6^ m2 defined 



3. Finally, the functions B generated by summing over all MHV Feynman diagrams 
with the range dictated by ([4.3)1 are precisely those included in the double sum for 
the finite box functions in the first line of ()2.2)1 (or (|2.19jl ) upon identifying m\ and 
m 2 with s and m. To be precise, (J4.3[) includes the case where the indices s and/or 
m (in the notation of ()2.2|) and (|2.19jl ) are equal to either i or j; but for any of these 
choices, it is immediate to check that the corresponding coefficient b l ^ s vanishes. 

This settles the agreement between the result of our computation with MHV vertices 
and ([2.19J1 for the part corresponding to the box functions. Next we have to collect the 
cuts contributing to particular triangles, and show that the z-integration reproduces the 
expected triangle functions from ([2.190 . each with the correct coefficient. 



in 023). 
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To this end, we notice that for each fixed triangle function T(p, P,Q), exactly four 
phase space integrals appear, two for each of the two possible cuts of the function. More- 
over, a gauge invariance similar to that of the box functions also exists for triangle cuts 
(see Appendix C), so that we can choose r\ in a way that the T numerators in 1)4.12)1 
become independent of z. A particularly convenient choice is 77 = hi, since it can be 
kept fixed for all possible cuts. Choosing this gauge, we see that a sum, T, of terms 
proportional to cut-triangles is generated from (|4.11|) (up to a common normalisation): 



T := T A + T B + T c + T, 



- D 



(4.15) 



where 
T c 



S(i,j,mi,m 2 ) _ S(i,j,m 1 - l,m 2 ) 
(mi • m 2 ) ((mi - 1) • m 2 ) 

S(i,j, m 2 ,mi) _ S(i,j,m 2 + I, mi) 
(m x ■ m 2 ) ((m 2 + 1) • mi) 

S(i, j,m 2 + I, mi - 1) S(i,j,m 2 ,mi-1 



{{m 2 + 1) ■ (mi - 1)) (m 2 • (mi - 1)) 
S(i, j,mi - l,m 2 + 1) <S(z, j,m 1 ,Tn 2 + 1 



((mi - 1) • (m 2 + 1)) (mi • (m 2 + 1)) 

Here we have defined 



S(i,j,m 2 ,P L )A A , (4.16) 
S(i,j,m u P L ) A B , 

S(i,j,mx - 1,P L )A C , 
S(i,j,m 2 + 1,P L ) A D . 



and A/, / 
A A 
A B 

A c 
A D 



S(a,b,c,d) = ta + (#ttJMMfci) , 
A, . . . , D, are the following cut-triangles, all in the PL ;z -cut 

1 



(4.17) 



(m 2 ■ P L - Z ) 
1 

(mi ■ P L]Z ) 
1 



((mi - 1) • P L , 
1 



((m 2 + 1) ■ P L . 



Q 2 -cut of - T(m 2 , P L;2 - m 2 , -P L -,z) , 
P 2 -cut of - T(mi, -P L -z, Pl,z ~ m x ) , 

- = Q 2 -cut of T(mi - 1, -Pl ;2 - (mi - 1), P L;2 ) , 

- = P 2 -cut of T(m 2 + 1, P L;z , -P L , Z - (m 2 + 1)) . 



(4.18) 



Next, we notice that the prefactors multiplying A B , Ac become the same, up to a 
minus sign, upon shifting mi — 1 — > mi in the second prefactor; and so do the prefactors 
of A a, A D upon shifting m 2 — > m 2 + 1. Doing this, — A B and the shifted A c become the 
two cuts of the same triangle function T(mi, —Pl jZ ,Pl-,z — m i), and similarly, — and 
Ad give the two cuts of the function T(m 2 , Pl- z — m 2 , —Pl- z )- Furthermore, in Appendix 
C we will show that summing the two dispersion integrals of the two different cuts of a 
triangle indeed generates the triangle function - in fact this procedure gives a novel way 
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Figure 5: A triangle function with massive legs labelled by P and Q, and massless leg p. 
This function is reconstructed by summing two dispersion integrals, corresponding to the 
P 2 - and Q 2 z -cut. 



of obtaining the triangle functions. 11 Specifically, the result derived in Appendix C is 



dz 

z 



{Pi 



+ 



(Ql 



(PzP) (QzP) 



2[7recsc(7re)] T e (p,P,Q) 



(4.19) 



where the e-dependent triangle function T e (p, P, Q) (with p + P + Q = 0) was introduced 
in ()2.14|) and gives, as e — > 0, the triangle function (|2.1Hj) (as well as the bubbles when 
either P 2 or Q 2 vanish). The result 1)4.19)1 holds for a generic choice of the reference 
vector 7), see (|C.7|) - (IC.12|) . We give a pictorial representation of the non-degenerate and 
degenerate triangle functions in Figures 5 and 6, respectively. 



At this point, it should be noticed that for a gauge choice different from r\ — k{ adopted 
so far, the numerators T in (|4.12j) do acquire an ^-dependence. This gauge dependence 
should not be present in the final result for the scattering amplitude. Indeed, it is easy 
to check that, thanks to (jC7|) . the coefficient of the ^-dependent terms actually vanishes. 

11 A remark is in order here. In our procedure the momentum appearing in each of the possible cuts 
is always shifted by an amount proportional to Z7j; the triangle is then reproduced by performing the 
appropriate dispersion integrals. Because of the above mentioned shift, we produce a non-vanishing cut 
(with shifted momentum) even when the cut includes only one external (massless) leg, say k, as the 
momentum flowing in the cut is effectively k z = k — zrj, so that k\ ^ 0. 
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Figure 6: A degenerate triangle function. Here the leg labelled by P is still massive, but 
that labelled by Q becomes massless. This function is also reconstructed by summing over 
two dispersion integrals, corresponding to the P 2 2 - and Q\-cut. 



Using now (|4.15j) - ()4.19j) and collecting terms as specified above, we see that the generic 
term produced by this procedure takes the form 



S(i,j,a,p m ) _ S(i,j,a + l,p m ) 

(k a ■ p m ) {k-a+l ' Pm) 

with P = <? a +i,m-i and Q = q m +i, a - 



S(i,j,Pm,Q) T(p m ,P,Q) 



(4.20) 



Finally, we implement the symmetrization of the indices i, j, as explained earlier, and 
convert (j4.20j) into 

c% a T(p m ,P,Q) , (4.21) 

where the coefficient c l £ „ is 12 



S(i,j,a + l,p m ) S(i,j,a,p m ) 



[ka+l ' Pr. 



(ka ' Pr, 



3i Vrai Qm,a) ~ 3i Qm,ai Pro) 



[(fcj + kj / 



212 



(4.22) 

which coincides with the definition of c l r ^ a given in ()2.6|) . Finally, it is easy to see that 
in summing over the range given by ()4.3|) . we produce exactly all the triangle functions 

12 In writing (|4.22Jl . we make also use of the fact that S(i, j,q m -i ia ,Pm) — S(i,j,q mta ,p m ). 
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appearing in the second line of 1)2.19)1 . It is also important to notice that the bubbles, 
which appear in the last line of (J2.2j) . are actually obtained as particular cases of triangle 
functions where one of the massive legs becomes massless, as observed at the end of 
Section 2. 

In conclusion, we have seen that all the terms in ()2.19j) . i.e. finite box contributions 
and triangle contributions - which include the bubbles as special (degenerate) cases - are 
precisely reproduced in our diagrammatic approach. 
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Appendix A: Passarino-Veltman reduction 



In Section 4 we saw that a typical term in the Af = 1 amplitude is the dispersion integral 
of the following phase space integral: 

C(m u m 2 ) := /dLIPSfe -Z x ; P L , Z ) tr 4^^Ji)tr+p^^ m2 r 2 ) (AJ) 
y (z • jj^mi • /i)(m 2 • / 2 ) 

The full amplitude is then obtained by adding the dispersion integrals of three more terms 
similar to (|A.1J) but with mi replaced by mi — 1 and/or m 2 replaced by m 2 + 1. The goal 
of this Appendix is to perform the Passarino-Veltman reduction [33] of (jA.ljl . which will 
lead us to re-express C(m 1 ,m 2 ) in terms of cut-boxes, cut-triangles and cut-bubbles. 

The explicit forms for the Dirac traces involve Lorentz contractions over the various 
momenta, so in a short-hand notation we can write these as 

T(z,j, mi )^^ := tr + (^ mi f x ). (A.2) 

C(mx, m 2 ) can then be recast as 

C(mi , m2) = T(i,i,m) r ni,i,m). inmumi]P ^ _ (A . 3) 

\ l ' 3) 

where 13 

X^(m 1 ,m 2 ,P i ) = [dUPS(l 2 ,-h;P L )- (A4) 

X M ^(mi, m 2 , Pl) contains three independent momenta mi, m 2 and P^. On general grounds, 
we can therefore decompose it as 

= TfJo + m^Ji + m^m v 2 X 2 + P£P£X 3 + m^m^X 4 

+ J 5 + r<P£Z 6 + P^lr + m^P[l 8 + P^m u 2 l 9 , (A.5) 

for some coefficients Ti, % = 0, ...9. One can then contract with different combinations of 
the independent momenta in order to solve for the Xj. For instance, two of the integrals 
that we will end up having to do are rj^I^ and m^m^T^. Using momentum conservation 
h~h J rPL — and the identity a-b = (a + b) 2 /2 = —(a — b) 2 /2 for a, b massless momenta, 
we can convert these integrals into ones which have the general form 

f(aM f dUPS(l 2 ,-h;P L ) 

J (h ■ mi y{l 2 ■ m 2 y ' lA -° j 

possibly with a kinematical- invariant coefficient, and with a and b ranging over the values 
1,0, —1. The results of these integrals are collected in Appendix B. As an example, we 
find that 

m«V = / dUPS(l 2 , -h; P L ) - K ■ P ^ I dUF f 2 >- ll > PL) . (A.7) 
J {h ■ m 2 ) J (l 2 ■ P L ) 



13 



For the rest of this Appendix we drop the subscript z in Pl- z for the sake of brevity. 
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Considering the values (a, 6), the case (1, 1) is a cut scalar box, (1,0) and (0, 1) are cut 
scalar triangles, (1,-1) and (—1,1) are cut vector triangles, whilst (0,0) is a cut scalar 
bubble. 



Because of the structure of T(i,j, and T(i, j, m 2 )j,, terms with coefficients such 
as T(i,j,mi)^T( y i,j,m 2 ) u mi tl m 2 u are zero, and thus some of the X, do not contribute to 
the final answer. The only contributing terms are found to be X 3 , X 5 , X 7 and X 8 , and we 
find that 



C(mi,m 2 ) 



X, 



+ 
+ 



tr+(^j^j^ mi -pi)tr_|_(^j|^j^ m2 ^ 

h)tr + 
(i-j) 2 

tr+(^i^j^ mi ^ m2 )tr_)_(^i^j^. 



X fi 



-*-7 



m 2 f L 



x s 



(A. 



The inversion of (|A.5|) in order to find the coefficients is tedious and somewhat lengthy, 
so we just present the results for the relevant Xj in (|A.8|) above: 



X 3 = ^{2(m 1 -m 2 )P L 2 X(°- )-A(m 1 -P i )X( 1 '°)+A(m 2 -P i )X( ' 1 ) 
+ 2(m 2 -P L ) 2 X(- 1 ' 1 )+2(m 1 -P L ) 2 X( 1 - 1 )} , 



(A.9) 



X, 



(mi • m 2 ) 2 N 2 
+ 3(m 1 -m 2 ) 2 (P|) 



4(mx ■ P L ) 2 (m 2 ■ P L ) 2 - 6(mx ■ P L )(m 2 ■ P L )(mi • m 2 )P 2 



J(0,0) + 



2(mi ■ P L ) 2 (m 2 ■ P L ) - -(mi ■ m 2 )P 2 



A(mi -P L )X (1 ' 0) 



2(mi • P L ) 2 (m 2 • P x ) - -(m a • m 2 )P| 



+ 2 



+ 2 



(mi • m 2 )P| - (mi • P L ){m 2 ■ Pl) 



(mi • m 2 )P| - (mi ■ P L ){m 2 ■ P L ) 



(m 2 -P i ) 2 X(-W 
(mi-Pi) 2 ^ 1 - 1 ) 



(A.10) 



X 7 



(mi • P/^) (mi • m 2 )N 2 



2 (mi • P L ) 2 (m 2 • P L ) 2 - 3(mi ■ P L )(m 2 • F z )(rm • m 2 )P 2 



• i m + \{mi- m 2 )PlN{m l ■ P L )X^ - (mi • P L )N(m 2 ■ P L fl^ 



2{ mi ■ P L )(m 2 • P L ) 3 X^ - (mi • m 2 )P 2 L {m x ■ P L ) 2 X^ 



-i) 



(A.ll) 
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x« 



(m 2 • Pl)(^i • m 2 )N 2 



2{m 1 ■ P L f{m 2 ■ P L f - 3{ mi ■ P L ){m 2 ■ P L )( mi ■ m 2 )P 2 L 



. j(o,o) + (m2 . p L ) N ( mx . p L fl^) - ( mi . m 2 )P|iV(m 2 • P L ) X^ 1 ) 
- (mi ■ m 2 )P!(m 2 • P L ) 2 X^ - 2{m x ■ P L f(m 2 ■ P L )l^' r ' 



(A.12) 



where N = (mi ■ m 2 )P 2 — 2(m\ ■ Pl)(t«2 ■ Pi). The explicit expressions for the relevant 
are summarised in Appendix B. 



Combining (jA.8|) and ()A.9|) - ()A.12|) with the identity (jD.7|) and the explicit expressions 
for the integrals X^ a '^ in Appendix B, we arrive at the final result (j4.12|) . 



Appendix B: Box and triangle discontinuities from 
phase-space integrals 



The integrals that arise in the Passarino-Veltman reduction in Appendix A have the 
general form: 

J {h-mxY{h-m 2 f ' 1 • ) 

where we have introduced dimensional regularisation in dimension D = 4 — 2e [34] in 
order to deal with infrared divergences. 

There are six cases to deal with: p lfi \ jC 1 ' 1 ), X^ 1 - 1 ), J (1 ' _1) , though due 

to symmetry we can transform X- 1,0 ) into X^ 0,1 ), and X -1 ' 1 ) into X^ x ~ x \ so we only need 
consider four cases overall. 

Generically we will evaluate these integrals in convenient special frames following 
Appendix B of [17], with a convenient choice for mi and m 2 . For instance, in the case 
of jC 1 ' 1 ) it is convenient to transform to the centre of mass frame of the vector 1% — l 2 , so 
that 

h = ^Pl ;2 (1,v) , l 2 = Ip L;2 (-l, v ), (B.2) 

and write 

v = (sin#i cos# 2 , . . . , cos#i) . (B.3) 
Using a further spatial rotation we write 

mi = {mi, 0, 0, mi) , m 2 = (A, B, 0, C) , (B.4) 

with the mass-shell condition A 2 = B 2 + C 2 . 
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After integrating over all angular coordinates except 6\ and 9 2 , the two-body phase 
space measure in 4 — 2e dimensions becomes [17] 



c^LIPS(Z 2 ,-Z i; P L;2 



i 

7T2" 



4r(|- e ) 



-2c 



d6 1 d6 2 (sinfli^-^sinfla) • (B.5) 



As a result of this and of our parametrizations of h, l 2 , ra% and m 2 , the integrals take the 
form 



7T2 



4 r(i- e ) 



L:z 



-2c 



J 



(0,6) 



where 



A (o,o) 
A (i,o) 

A (o,D 

A (i,D 

A (-i,i) 
A (i,-D 



;b.6) 



;bj) 



and J^ a ^ is the angular integral 

«2tt 



J 



(0,6) 



2 
4 

-mi , 
-"^2 , 



(sin^O^fsin^)-^ 



1 J 2 {l-cos6 1 ) a {A + Ccos6 l + Bsm6 1 cos6 2 ) 



;b.s) 



The integrals (jB.8|) have been evaluated in [35] for the values of a and b specified above, 
and we borrow the results in a form from [36]: 



j(0,0) 
j(-lA) 



2tt 



1 - 2e 

7T 

■ 



(B.9) 



e 

7T 1 / A-C 

— - 2J P! 1,1,1 -e, 



e A 
2tt(1 - e) 
efl -2e) 



2A 



2 Pi 



A-C 



-1,1,1-e, 

' 2A 



Here, A and P will differ depending on which case we are considering and our particular 
parametrization for it, but in all cases the combinations that arise can be re-expressed 
in terms of Lorentz-invariant quantities using suitable identities. In the case of JP 1,1 ) for 
example, one uses the easily verified identities 



N(P L , Z ) = -Pl z (A + C) mi , m 1 -m 2 = m 1 (A-C) 



(B.10) 
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where N(Pl ;z ) was defined in (J4.14I 



Eventually, after re-expressing A and B in this way, and upon application of some 
standard hypergeometric identities we find the following: 



A -lj(0,0) 
A -lj(l,0) 

A -ij(o,D 

A -1 J (1,1) 
A -ij(-i,i) 



2ty 



l-2e ' 
1 2tt 



e rri! ■ P L]Z 
1 2tt 
e m 2 • P L - Z ' 

Svr /I / (mi ■ m 2 )P| ;2 



7T 



(mi • P L . Z ) 



N(P, 



L:z, 



+ t( mi ' p i;-)( m 2 ■ ^l ;2 ) - (mi ■ m 2 )Pf J J 



A -ij(i,-D 



7T 



K • p L , z y 

where A is the ubiquitous factor 



7T2 



P, 



-2e 



4r(i-e) 



(B.ll) 



(B.12) 



Appendix C: Reconstructing triangles from 
dispersion integrals in a gauge- invariant way 



In this Appendix we find a new representation of the triangle function 

T{ ^ P , Q) = ^m, ( c.D 

as the dispersion integral of a sum of two cut-triangles. 14 A comment on gauge- (in) dependence 
is in order here. Recall from Section 2, Eq. ()3.1|) . that in the approach of [17] to loop 
diagrams one introduces an arbitrary null vector 77 in order to perform loop integrations. 
The corresponding gauge dependence should disappear in the expression for scattering 
amplitudes. In what follows we will work in an arbitrary gauge, and show analytically 

14 For a review of dispersion relations, see [37]. 
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that gauge-dependent terms disappear in the final result for the triangle function. Perhaps 
unsurprisingly, this gauge invariance will also hold for the finite-e version of T(p, P,Q), 
which we define in ([2. 14)1 . 

To begin with, recall from 1)4.18)1 that the basic quantity we have to compute reads 



K 



dz 



(PzP) (QzP) 



where P + Q + p = 0. We will work in an arbitrary gauge, where 

P z := P - zt) , Q z := Q +zr] . 
A short calculation shows that 

P z p = Pp 
QzP = Qp 

where 



i - b P (p 2 - Pi) 

l-b Q (Q 2 -Q 2 z ) 



r]p 



VP 



' ■ 2( V P)(pP) ' 
It is also useful to notice the relation 



2(vQ)(pQ) ' 



On bp 



(C.2) 

(C.3) 

(C.4) 
(C.5) 

(C.6) 
(C.7) 



as well as (Pp) = —{Qp) = (1/2) (Q 2 — P 2 ), which trivially follows from momentum 
conservation. We can then rewrite (IC.2I) as 



where 



1 



(Pp) 
7r cscfVe 



(Pp) 



n = ii - x 

ds'(s')- e 

-P 2 )-' 



2 ■ 



1 



s' - P 2 )[l-b P (P 2 - s')] 

-bp 



(C.8) 



(C.9) 



b P P 2 - 1 



Xo 



But ()('.7)1 implies 



1 



(PP) 
7T CSC(7T€ 
~JPpJ 



ds'is 



1 



-Q 



2\-e 



s> - Q 2 )[l - b Q (Q 2 - s')] 
-bo 



b Q Q 2 - 1 



b P P 2 - 1 b Q Q 2 - 1 
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(CIO) 



(C.11) 



so that we can finally recast ()C.2j) as: 

K = 2[ 7 recsc( 7 re)]l ( ~ jP2) g2 ~^; g2) ^ = 2 [ne csc(vr e )] T e (p, P, Q) , (C.12) 
where the e-dependent triangle function is 15 

Te(p,P,Q) ■= ^X^p^ ■ (C13) 

This is the result we were after. Notice that all the gauge dependence, i.e. any dependence 
on the arbitrary null vector 77, has completely cancelled out in (|C.12[) . 

We now discuss the e — > limit of the final expression (jC.12j) . As already discussed 
in Section 2 (see (j2.15|) and (|2.16jl ). in studying the e — > limit of TZ (and hence of 
T e (p, P, Q)) we need to distinguish the case where P 2 and Q 2 are both nonvanishing from 
the case where one of the two, say Q 2 , vanishes. In the former case, we get precisely the 
triangle function T(p, P, Q) defined in (jC.lj) : 

YmxK = 2 T(p,P,Q), P 2 ^ , Q V . (C.14) 
In the latter case, where Q 2 = 0, we have instead 

hmTZ = - 2 - { ~ P p} - , PV0,Q 2 = 0, (C.15) 

which corresponds to a degenerate triangle. 

The final issue is that of the gauge invariance of the contributions to the amplitude 
from the box functions B (this is also relevant to the issue of gauge invariance in the 
Af = 4 calculation of [17], and in that paper we also gave a general argument for gauge 
invariance). We expect that an explicit analytic proof of the gauge invariance of the box 
function contribution to the amplitude could be constructed using identities such as those 
in Appendix B of [17]. In the meantime, numerical tests have shown that gauge invariance 
is present [30]. Indeed, it would be surprising if this were not the case, given that the 
correct, gauge invariant, amplitudes are derived with the choices of gauge we have made 
here and in [17]. We have also carried out the MHV diagram analysis of this paper using 
the alternative gauge choice r] = k m2 ; one obtains (|2.19|) . 



Appendix D: Spinor and Dirac-trace identities 

In this appendix we present some identities that are needed for many of the manipulations 
required to massage the M = 1 amplitude into the form given by BDDK 12. 2[) . 

15 The e-dependent triangle function already appeared in H2.14|l . 
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For spinor manipulations, the Schouten identity is very useful: 

(ij)(kl) = (ik)(jl) + (il)(kj) . (D.l) 

Furthermore we have: 

(tj)[jt} = ti + (Mi) = 2(fci • kj) , (D.2) 
(i]){]l}{lm)[mi] = ti + {UMm) , (D.3) 

(ij)[jl](lm)[mn](np)[pi] = tT: + {$$$$ m $ n $ p ) , (D.4) 

for momenta fcj, fcj, fcj, fc m , k n , k p . For a nice introduction to the spinor helicity formalism, 
see [38]. 

For dealing with Dirac traces, the following identities are useful: 

to+(ftiftj$iftm) = 4(/cj • kj)[ki ■ k m ) — tT+(pjpifiip m ) , (-D-6) 

for similarly generic momenta. If fcj, kj, k mi and fc m2 are massless, while Pl is not neces- 
sarily so, then we have the remarkable identity: 

2(/c mi • i 3 i)tr+(^j^j^ mi ^ m2 )tr+(^j^j^ m2 ^'^) 
—2(k m2 • i\)tr+(^j^j^ mi ^£)tr+(^j^j^ m2 |c mi ) = . (-D.7) 
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